Anomalous specific heat jump in a two-component ultracold Fermi gas 
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The thermodynamic functions of a Fermi gas with spin population imbalance are studied in the 
temperature-asymmetry plane in the BCS limit. The low temperature domain is characterized by 
anomalous enhancement of the entropy and the specific heat above their values in the unpaired 
state, decrease of the gap and eventual unpairing phase transition as the temperature is lowered. 
The unpairing phase transition induces a second jump in the specific heat, which can be measured 
in calorimetric experiments. While the superfluid is unstable against a supercurrent carrying state, 
it may sustain a metastable state if cooled adiabatically down from the stable high-temperature 
domain. In the latter domain the temperature dependence of the gap and related functions is 
analogous to the predictions of the BCS theory. 



Recent experiments P, on ultracold dilute gases of 
fermionic atoms trapped unequal number of fermions in 
two different hyperfine states. These experiments started 
addressing some of the long standing problems in the 
theory of asymmetric superconductors (ASC) that are of 
interest in variety of fields includingrnetallic supercon- 
ductors LJ. Lul . nuclear systems [a, H 13 and high density 
QCD [gl Tl oTfl lL fl^| . The unprecedented control over the 
many-body systems achieved in the experiments with ul- 
tracold dilute fermions combined with the possibility of 
tuning the interactions via the Feshbach resonance mech- 
anism provide for the first time a realistic perspective 
of testing the predictions of the theories of ASC in the 
context of dilute fermionic systems. The realizations of 
various phases of ASC of dilute fermions have been in- 
tensively studied on the theoretical front; the simplest 
realizations are the isotropic, homogeneous phases that 
are characterized either by a Zeeman splitting of Fermi- 
levels ITS Il4l |l5l or by pairing between light and heavy 
fermions |lq . Il7| . At large asymmetries the phases with 
broken space symmetries j^lTll2fjll2ll2ll23ll2ll25t 
and the mixed phases |2(| [2g become energetically more 
favorable. Alternatives include pairing in higher angular 
momentum states 28, 29] . Finite size and trap geometry 
introduce an additional complication to the problem and 
may qualitatively affect the comparison between the the- 
ory and experiment [3(], El ■ A number of related prob- 
lems of interest are the nature of phase transitions be- 
tween the various phases and their relation to the topol- 
ogy of Fermi-surfaces 132. 1331 as well as the features of 
the BCS-BEC crossover |34U35ll36| j under population im- 
balance. 

The population asymmetry in ASC can be character- 
ized either in terms of the difference (mismatch) in the 
chemical potentials or the difference in the densities of 
the species. The first case arises when the "chemical" 
equilibrium between populations admits transmutation 
between the different spin states, as e. g. under the 
equilibrium with respect to the weak interactions in cold 
dense hadronic/quark matter. We shall specify our dis- 



cussion from the outset to the second case and assume 
that the total number of fermions is fixed and the asym- 
metry is maintained with respect to the number densities 
of different species, as is the case in the experiments on 
ultracold fermions. 

The aim of this Letter is the study of the temperature- 
asymmetry phase diagram of an ultracold Fermi gas 
with pairing between two unequally populated hyper- 
fine states in the BCS regime. We propose that calori- 
metric experiments, which are within the current ex- 
perimental capabilities [33], can test the specific fea- 
tures of the phase diagram discussed below. The crit- 
ical temperature of metallic superconductors, according 
to the Bardeen-Cooper-Schrieffer (BCS) theory, is given 
by T c = 1.14 wrj e~ 1 / v , where lod is the Debye frequency 
and V is the dimensionless interaction. The parameters 
in this equation can be determined from a single calori- 
metric experiment; T c is determined from the position 
of the jump in the specific heat, and ljd from the slope 
of the specific heat over T 3 in the limit T — ► 0. The 
case of ASC is complicated by the fact that the asymme- 
try leads to a loss of coherence at low temperatures and 
the critical temperature and pairing gap become com- 
plicated functions of spin imbalance. This paper studies 
the impact of the asymmetry induced decoherence on the 
thermodynamics of ultracold gases; in particular we show 
that two jumps in the specific heat of ASC are possible 
if the asymmetry is large, albeit the second anomalous 
jump occurs withing the temperature domain where the 
superfluid is in a metastable state. The anomalous jump 
in the specific heat is a manifestation of the reentrance 
effect, i. e., the restoration of pair-correlations at finite 
temperatures T < T c 0, |3j| |33 • Below, we shall con- 
fine ourselves to the case of infinitely extended systems. 
Since the experiments are carried out in finite geome- 
tries, finite-size corrections need to be taken into account 
in a more complete analysis. Recent experiments 0, Q 
which have measured the density profiles of trapped gases 
with population imbalance demonstrate the importance 
of these effects in determining the ground state structure 



of the condensate |3fJ, |31J . Related work on finite tem- 
perature phase diagram of ASC appeared in refs. {IclE!) 
while our work was in preparation/revision. 

We approximate the pairing interaction by a zero-range 
force, which is characterized by the s-wave scattering 
length as- Our discussion is specified to the case where 
two hyperfine states of 6 Li are populated (the scattering 
length in units of the Bohr radius is as/aB = —2160, 
but can be varied at will via Feshbach resonance mecha- 
nism). The fermion masses are assumed to be equal; the 
extension to the case of unequal masses (as would be the 
case in the mixtures of 6 Li and 40 K) is straightforward. 

We consider a uniform gas of fermionic atoms in two 
hyperfine states (spins) labeled as j and |; the interaction 
Hamiltonian is 



Hint. = -V 



a/3 



where V>£(a:) and i> a (x) are the creation and annihila- 
tion operators of a state at the space point specified 
by the position vector x and spin a(=1,l) and V is 
the two-body bare contact interaction. The normal and 
anomalous propagators are [2(j Gn(p) = oj + £s±£a/D, 
F\p) = At/D , where D = (ui — £ A f - £| - A 2 and 
£s = ( e pT + e p|)/2 an d £a = (£p| — e pl)/2 are > respec- 
tively, the parts of the quasiparticle spectrum which are 
symmetric and antisymmetric under time-reversal oper- 
ation, e p f , e p i , are the single particle energies in states 
| and [. The dispersion relation of the quasiparticles in 
the paired state is w-m = £a ± ^/£§ + A 2 . The limit 
£a — * corresponds to the case of equal number of spin- 
up and down particles. The explicit form of the symmet- 
ric and antisymmetric (under time-reversal) parts of the 
quasiparticle spectrum are £s = (p 2 /m — p^ — /x*)/2 and 
£a = (— jU|+/iT)/2, where m is the atom's bare mass, p is 
the relative momenta of fermions bound in a Cooper pair 
in the center-of-mass frame at rest, the effective chemical 
potentials (J**n = A 1 !/! — ^i/T include the constant shift 
due to the self-energy Sf/j. = Topm, where To is the 
two-body scattering T-matrix. The gap equation is 
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where Uo is the strength of a contact interaction and A 
is an ultraviolet cut-off (for details see Hjj). The occu- 
pation probabilities of species are given by 

n m (p) = u% [f(u m ) - /(-w i/T )] + f(-u> in ), (2) 



where u 2 p = \/2+£ s /{2^/£ 2 s + A 2 ); these are normalized 
to the densities of species pyi = J2 P n Vl(p)- The free- 
energy is T = E P ( £ pT n p,T + e Pl n P,l) ~ u o lA2 - TS s- 
The free-energy of the normal state follows by setting 
in this expressions A = = Se. The entropy of ASC 
is defined in terms of the temperature derivative of the 
free-energy Ss — —dTjdT. The specific heat follows as 
C v = T(dS s /dT) = -T{d 2 T/dT 2 ). 




FIG. 1: Dependence of the pairing gap (upper panel) the free- 
energy difference (middle panel) and the superfluid density 
(lower panel) on the temperature for several values of the 
density asymmetry. The instability domain (p s < 0) lies to 
the left of triangles. 



The local stability requires that the free energy is a 
convex function of the appropriate variables and it has 
been established that homogeneous ASC could become 
unstable in this sense 0, l25l l2rj 

sum in mm 

Specifically, (A) the system is unstable against phase sep- 
aration unless the curvature matrix Xij — 
is positive definite. This implies that either the eigen- 
values Ai, a = Trxv/2 ± VC 1 ^) 2 ™ 4DetXy/2 > or, 
equivalently, xtt > 0> Xij > (Sylvester criterion). 
Further, (B) the condition d 2 Ts/dA 2 > needs to be 
fulfilled. Finally, (C) the system may become unstable 
against spontaneous generation of currents (formation of 
the LOFF phase) when \q — d 2 J r s/dq 2 < 0, where q is 
the center-of-mass of momentum of a Cooper pair. The 
latter instability manifests itself in the negative super- 
fluid density, p s = m 2 x{q)\ q =o 1 and purely imaginary 
Meissner mass. For small quasiparticle velocity u the 
leading order contribution to the ratio of the superfluid 
to the total density is [44| 
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We fix the density of 6 Li atoms at p 



3.8 x 10 



12 



cm , which in the case of a spatially uniform and spin- 
symmetric system at T = corresponds to the Fermi 




FIG. 2: Entropy (upper panel) and the heat capacity (lower 
panel) as a function of the temperature for several values of 
the density asymmetry. Note that for a — 0.02, cv(T = 0) = 
4.5 is finite. The grey lines show the values in the normal 
state. 



momentum fej? = 4.83 x 10 4 cm -1 and kpas = —0.558. 
Since we work in the BCS regime the conditions (A) and 
(B) are fulfilled (in agreement with ref. [iof). Figure 1 
displays the A(T), AJF(T) = F N -F S and p s (T)/p func- 
tions for constant asymmetry a = (p-\ — Pi)/(p-\ + Pi)- 
The length and energy are measured below in units of 
L = 1CP 4 cm and e = 10 2 nK, unless explicitly spec- 
ified. Near the critical temperature the asymptotic be- 
havior of the pairing gap for T — > T c is described, to lead- 
ing order in (A/T) 2 , by the BCS-type relation A(a) ~ 
[T c (a)(T c (a) — T)] 1 / 2 , i. e., the high temperature portions 
of the A(T) curves are self-similar. The low-temperature 
BCS asymptotics A(T) - A(0) - v /2 7rTA (0)e" A(0)/T is 
qualitatively modified even at small asymmetries, since 
the gap decreases as T — > 0, instead of staying constant. 
At large asymmetries (e.g. a > 0.04) the reentrancc 
effect sets in: the gap is non-zero only in a finite do- 
main of temperatures bounded by two critical temper- 
atures. The physical origin of the upper critical tem- 
perature is analogous to the BCS case, where the pair- 
ing correlations are destroyed by thermal fluctuations. 
The lower critical temperature is due to the loss of co- 
herence induced by the asymmetric population and is 
specific to ASC. At large asymmetries the temperature 
dependence of the gap is A (a) ~ [T*(a)(T - T^a))] 1 / 2 , 
where T* is the lower critical temperature. The free- 
energy plots (Fig.nj reflect the temperature dependence 
of the condensation energy, which scales as A 2 (a, T); the 
temperature dependence and asymptotics of AJ-(T) is 



FIG. 3: Occupation probabilities of the majority and minority 
components as a function of the momentum for several asym- 
metries at temperatures T = 1 nK (upper panel) and T — 25 
nK (lower panel) . All lines correspond to stable states, except 
a 7^ lines in the upper panel. The grey line shows the same 
for the unpaired symmetric state at T — 0. 



understandable in terms of this scaling. In the temper- 
ature domain where p s {T) < 0, the homogeneous ASC 
is metastable; there exists a lower extremum (perhaps 
minimum) of the free energy corresponding to the cur- 
rent carrying LOFF phase. Nevertheless if prepared at 
high enough (but T < T c ) temperature and cooled down 
adiabatically, the metastable phase can be sustained long 
enough to carry out measurements. The temperature de- 
pendence of the entropy and specific heat (more precisely 
C v /T) is shown in Fig. [5] At T — > T~ the entropy scales 
linearly with temperature, Ss oc T — T c , with asymme- 
try dependent slope. In the low-temperature metastable 
region the superfluid entropy is anomalous, since its ab- 
solute magnitude is larger than the entropy of the normal 
state, i. e., the superfluid appears to be less ordered than 
the unpaired state. The temperature for the onset of 
anomalous regime (Ss > Sn) coincides with that for the 
onset of instability within our numerical accuracy. 

The ratio AC/T C , where AC = Cs — Cn is the jump in 
the specific heat at the critical temperature, depends only 
on the density of states and is a universal characteristic 
of a system. The jump itself is a characteristic feature of 
a second order phase transition that allows to locate T c 
experimentally. There is a second (anomalous) jump in 
the specific heat associated with the reentrance effect at 
the lower critical temperature T*, which lies within the 
metastable domain. At T — > T~ the specific heat scales 
as Cv oc \T— T c \. Its low-temperature asymptotics differs 
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dramatically from the predictions of BCS theory, where 
the superfluid "thermal inertia" is small compared to the 
normal case. The specific heat of ASC is larger than 
that of the normal state at sufficiently low temperature 
for any asymmetry. Calorimctric experiments aimed at 
measuring the specific heat of ASC can locate the critical 
temperature and observe the reentrance effect through 
the second (anomalous) jump ACy, if the system can be 
maintained long enough in the metastable state. 

The occupation probability of the majority and mi- 
nority components are shown in Fig. [3] for the high- and 
low-temperature regions of pairing. The remarkable dif- 
ference between these two cases arises from the fact that 
in the low-temperature regime the minority component 
is excluded from the region around the Fcrmi-momentum 
("blocking region" 6] or "breach" (lallll)- The deple- 
tion is large for large asymmetries. In contrast, the high- 



temperature regime does not feature a depletion region 
and the occupation are smooth functions of the momenta 
for arbitrary asymmetries. 

The finite range of interactions and finite size of the 
systems are not likely to modify the conclusions reached 
here. Indeed, reentrance behavior(s) have been predicted 
in paired nuclear systems which are characterized by 
complex finite range interactions 0, l38| and ultra-small 
metallic grains, which contain a small number of fermions 
and a single (odd) unpaired particle 39]. 
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